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Recently methods have been developed to extend the resummation of large-x double logarithms in 
inclusive deep-inelastic scattering (DIS) to terms not addressed by the soft-gluon exponentiation. 
Here we briefly outline our approach based on fixed-order results, the general large-jc structure in 
dimensional regularization and the all-order factorization of mass singularities, which is directly 
applicable also to semi-inclusive e+e^ annihilation (SIA). We then present some main results 
for the corresponding timelike splitting functions and transverse and longitudinal fragmentation 
functions. The close relation between DIS and SIA facilitates the determination of additional 
third-order results for the latter function which is fully known only at the next-to-leading order 
Therefore all above quantities can be resummed at next-to-next-to-leading logarithmic accuracy. 
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1. Introduction 

In the past years quite a few studies have addressed the threshold behaviour of higher-order split- 
ting functions and hard-process coefficient functions in perturbative QCD beyond the quantities or 
contributions covered by the standard soft-gluon exponentiation [1-13]. Except for the gluon con- 
tributions to the longitudinal structure and fragmentation functions Fj^ which include an additional 
factor (1— x), the dominant terms are of the divergent but integrable 'double-logarithmic' form 

a;'ln2"-"o-f(l_;c) . (1.1) 

Here depends on the quantity under consideration, with nQ = 2 for the 'off-diagonal' splitting 
functions Pij{x) with i j - recall that the diagonal splitting functions Pii{x) do not exhibit any 
n-dependent large-x enhancement in the MS scheme adopted in the present contribution [14, 15] - 
and nQ = l for the structure functions in deep-inelastic scattering (DIS) and fragmentation functions 
in semi-inclusive e^e^ annihilation (SIA) with the exception of the functions Fl where = 2. 
£ = 0, 1, 2, ... in Eq. ( |1.1[ ) represent the leading logarithmic (LL), next-to-leading logarithmic 
(NLL), next-to next-to-leading logarithmic (NNLL = N^LL), . . . large-x contributions. 

Using the third-order results of Refs. [16, 17], the LL, NLL and NNLL coefficients of the 
(non-singlet) quark coefficient functions for the structure functions F2, Ft, and Fl in DIS and their 
counterparts Fj, Fa and Fl in SIA have been obtained in Refs. [8,9], together with the correspond- 
ing LL and NLL results for the Drell-Yan process, from the single-logarithmic behaviour of the 
physical evolution kernels of these observables. The corresponding approach is not sufficient for 
all-order predictions in flavour-singlet cases, but was used to obtain the NNLL approximations of 
the fourth-order contributions to the off-diagonal spacelike splitting functions for the parton 
distribution and the longitudinal coefficient function Ci.g in DIS [10, 1 1]. 

Those DIS results have been confirmed and extended by studying the Z)-dependence of the 
unfactorized structure functions in dimensional regularization [12, 13]. This approach facilitates 
the determination of the one previously missing parameter in the N^'LL non-singlet coefficients 
for F2 and Ft, - a. result that has been obtained independently in Ref . [6] - and the extension of the 
NNLL all-order resummation to the off-diagonal splitting functions and gluon coefficient functions. 
In the present contribution we briefly report on the generalization of those results to the off-diagonal 
timelike splitting functions Pfj for parton fragmentation and the singlet coefficient functions for the 
fragmentation functions in SIA; a detailed account of our results will be published elsewhere [18]. 

For this purpose we consider the transverse and longitudinal gauge-boson exchange fragmen- 
tation functions Fj and as defined in Ref. [19] together with the corresponding Higgs-decay 
quantity in the heavy top-quark limit F^ . Schematically our notation for the corresponding time- 
like coefficient functions Cj, and the evolution of the fragmentation distributions at the physical 
scale with as = cCs{Q^) / {4n) reads (note the transposition of the splitting-function matrix) 

FjiQ^) = Cl{as)®DiQ^) , -j^Dj = P^j{as)®Di (1.2) 
where (g) represent the Mellin convolution and the appropriate summations over / are understood. 
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2. Outline of the resummation 

The primary objects of the present resummation are the unfactorized SIA fragmentation functions 

Plk = Caj(^Zl for a,k = T,g, ^,q, L,q and L,g (2.1) 

in D = 4 — 2£ dimensions. The functions Cj,- are given by Taylor series in e, with the terms 
including k more powers in ln(l— ;c) than the 4-dimensional coefficient functions. The timelike 
transition matrix consist of only negative powers of e and can be written in terms of 

P'^ = a,PQ + oIp^ +a^P2 + ... (2.2) 

and the corresponding coefficients j3„, of the beta function. This dependence can be summarized as 

Po,po , <£-"+^: +PuPi , ••• , : Pn-i ■ (2.3) 

Hence fixed-order knowledge at N'"LO (i.e., of the splitting functions to P^ and the corresponding 
coefficient functions) fixes the first m+l coefficients in the e expansion of FJj^ at all orders in as- 
The large-;c expansions of FJ_^j^ ^ (the corresponding relation for Fi is slightly different) are given by 



^^\a"e-'^''- = ^i,?ln"+^-i(l-x) + ^fyin"+^-2(i_^)^^_^ _ (2.4) 



If the constants up to ^j™' are known for all n and I, then the splitting functions and coefficient 
functions can be determined at N'"LL accuracy at all orders of the strong coupling. 

As in DIS, the order large-;c contributions to FJ^^ ^ are built up from n term of the form 

(A„,,£-2"+i+B„,,£-2«+2 + C„,,£-2"+3 + ...)(l-^)-*^^ , k=\,...,n (2.5) 

which arises from the phase-space integrations for the undetected final-state partons and the loop 
integrals of the virtual corrections [21-23]. Since the terms with £^^"+^, . . . , e^"^^ have to cancel 



in sum (2.1), there are «— 1 relations between the LL coefficients A^k which lead to the constants 



^jl^^ in Eq. ( |2.4| ), n — 2 relations between the NLL coefficients B„ k etc. As discussed above, a 



N'"LO calculation fixes the (non- vanishing) coefficients of £ ",...,£ at all orders n, adding 



m + l more relations between the coefficients in Eq. (2.5). Consequently the highest m+l double 
logarithms, i.e., the N'"LL approximation, can be determined in this manner from the N"'LO results. 

The coefficient functions for Ft = Fj , FJ^ and F^ are known at the second order [23-26]. 
The third-order timelike splitting functions have been determined, up to an uncertainty which is 
irrelevant in the present context, in Refs. [26-28]; see also Ref. [7] for their large-;c logarithms. 



The all-order factorization of the quantities (2.1) requires corresponding large-x results for the 



quantities Fj ^ and F^ ^. These are available from the soft-gluon exponentiation to a far higher 
accuracy than required here, cf. Ref. [29]. The calculations are carried out in Mellin-A^ space with 

\xv^{\-x) ^{-\)^N-\\xv^N ^\k(k-\)C,2\xv^-'^N ^ ...) , lnA^ = lnA^+7e (2.6) 

where /e is the Euler-Mascheroni constant. The required formalism is completely analogous to that 
in Ref. [13]. Our symbolic manipulations have been performed using FORM and TFORM [30, 31]. 
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3. Results for the timelike splitting functions 

As their spacelike counterparts, the LL and NLL contributions to the off-diagonal timelike splitting 
functions can be expressed in terms of functions S^n{^) introduced and discussed in Refs. [12, 13], 

^ tin „ . „ , v-l ti„ 



t(x) 



and 



(3.1) 



for = 0, 1 , 2, . . ., where S„ are the Bernoulli numbers in the normalization of Ref. [20]. We obtain 



(12Cf -6j8o) ^=^-i(-as) - ^^-2{-a,)+{6CF-po) 



+ NNLL contributions + ... , 
+ NNLL contributions + 



(3.2) 



(12Cf -2j8o) ^=^-i(as) + ^^-2{as) + {8Ca -ICp - Po)^i {a. 



(3.3) 



Here Ca and Cf are the standard SU (ric) colour factors with Ca = Uc = 3 and Cf =4/3 in QCD, 
rij- represents the number of effectively massless quark flavours, jSq = 11/3 Ca — 2/3 is the first 
coefficient of the (four-dimensional) beta function, and we have used the shorthand 



as = 4as Caf In with Caf =Ca-Cf 
reflecting the vanishing of the double logarithms for the 'supersymmetric' case Ca = Cf- 



(3.4) 



The first and second lines of Eqs. (3.2) and (B3b are the LL and NLL results, respectively. 



They differ from their spacelike counterparts, n^/C/r/'gq and Cf /?ijPqg only by coefficients of 
proportional to Ca — Cf- These differences are due to different (1— e)^^ pref actors of the spacelike 
and timelike unfactorized structure functions, e.g., the absence of the DIS gluon spin-averaging 
in the corresponding SIA quantity, which are otherwise related by a (at this level) simple analytic 
continuation. As in the spacelike case, we have no closed all-order expression for the NNLL 
contributions to Eqs. 



p and ( |3.3[ ) which therefore will be presented via tables to a sufficiently 
high order in ttj in Ref. [18]. After transformation back to ;c-space the fourth-order results read 



4w 



\n^{l-x)Cii,nf 



+ ln'^{l-x)CAFn 



22 ^ 
27 ^AF 

1432 
81 



27 

64 
9 



27"/ 



AF + 



1471 
54 



8C2) CafCf 



(3.5) 



- ^ CAFftf -^C^ + S CFUf + If + ^ (ln^(l -^)) 
= \n\l -x)C^pCf 
+ ln'^(l-:ic)CAFCf ( ^ — ^ (^2) C^F + ( ~ ^^2) CafCf 



26 ^ _ 14 ^ _ ±„ 

27 ^AF 77 77 



27 



27 "/ 



(3.6) 



212 r „ 

-gj- CAFnf 



13 ^2 
81 



17 



jCfUj 



81 "/ 



+ ^(ln^(l-x)) 



The LL coefficients vanish at all even orders in ttg from the fourth due to B2«+i = for n > 1. 
The above results are illustrated in Fig. |l| below for ttj = 0. 12 and rij^ = 5, i.e., at a scale ~ MJ. 
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Figure 1: The relative leading-logarithmic (LL), next-to-leading logarithmic (NLL) and next-to-next-to- 
leading logarithmic (NNLL) higher-order large-jc corrections to the third-order off-diagonal time like splitting 
functions Pj in Mellin-A^ space at a typical high-scale reference point. 
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Figure 2: The absolute LL, NLL and NNLL higher-order threshold corrections to the second-order trans- 
verse and longitudinal gluon coefficient functions Cr.g and in A^-space at a scale 2^ — M^. Note that our 

i 

2 

large-x limits are CT.g{x,as 



normalization of both functions differs by a factor of j from that of Refs. [23-25], i.e., here the lowest-order 

2Cf fls In (1 — x) + . . . and C/^,(ji:,fls) = ACpa^ (1 — x) + . . . . 
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4. Results for the SIA coefficient functions 

The NNLO results facilitate the resummation of also the (more complicated) coefficient functions 
Cr_g and C^^ at NNLL accuracy. Here we show, for brevity, only the NLL resummation of the 
former quantity (the NNLL contributions and the result for the latter will be given in Ref. [18]), 

1 Cf 



NCT,g{N,Ok 



2\nN Ca-CfI 
1 Cf{3Cf-^o 



fls Cf 
'4 Ca-Cf 



exp(2asCi7 In^A^)^o(as) - exp(2asCA In^A^) 



(4.1) 



exp (IusCa In^N) (8Ca + 4Cf - Po 



Cf 



3 Ca — Cf 

fls Cf 



Ca exp(2asCA In^A^) - Cf exp(2asCf hi^N)^o{as 
exp {IgsCf In^N) \ - GCF-^oia^) - (8Ca - 2Cf -po)^i (^s) 

+ NNLL contributions + 



4 Ca-Cf 

- (12Cf-4/3o)^^-i(as)- ?^^-2(a 
CIq etc 



Also this expression differs from its spacelike counterpart in Ref. [13] only in the coefficient of i 
The third-order contribution, now including the NNLL term, is given by 



CtAx) 



= \n^{\-x)CF 

+ ln'^(l-x)CF 

+ ln^(l-x)Cf 
41 112 



2^2 I 10^2 



1 n „ 269^2 , 17^ „ 338^ ^ 97^2 
jl^Aii-f ~ 54 t-A "T 21 F f ~ 27 ^F^^ ~ X^yp 

2990 16 ^ \ ^2 , / 3652 88 r \ ^ r 
-gl -g'sljCA + I -gi -^L,2)CfLa 



^-^C2)c|-^CA«^.-^C^n^. + ^(ln2(l-^)) 



436, 



(4.2) 



For the longitudinal fragmentation function the second-order results represent only the NLO 
contribution. The resulting NLL results for the gluon coefficient function read 

N^Cl^{N,a,) = 4a,CFexp(2a,CAln^N) + 2a,CFNC^^{N,a,) 



+ 8a; InNuf exp{2a,CA In'^N) ^ [4Ca - Cf j + ^ In-'A^CAjSo . 
where the first terms is the LL contribution. The resulting third-order x-space expression is 

= 8CFCiln^(l-x) 



il-x)-'clJx) 



+ ln3(l-x)CF 



20^2 I 52 ^ ^ 952 ^2 I 16^ „ 



+ 



(4.4) 



The third logs for Cl^ can not be derived by resumming the NLO results. The coefficient at order 
a^, however, can be obtained be comparing the physical kernels for (Ft, F^), cf. Ref. [32], to the 
analogous DIS results [11, 13] along the lines of Refs. [7,26], yielding the continuation of Eq. 



. + \n'{l-x)CF (62 - 32 C2) C^^ - ( ^ - 32 C2 ) CfCa + ^Cl + fCFnf-^ Ca 



5720^2 , 16, 



224 , 



+ ... 



where, as in Ref. [11], an additional d"^'''dahc contribution has been suppressed for brevity. The 
derivation of this result and its extension to the ln(l— x) term will be discussed in Ref. [18]. The 
numerical size of the resummed large-A^ corrections is shown in Fig. 2, again using ~ MJ. 
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5. Discussion and Outlook 

We have derived the all-order resummation of the three highest (NNL) threshold double logarithms 
for the off-diagonal timelike splitting functions and the coefficient functions Cf ^ , C/^^ , C'^ ^ and 
q gauge-boson and (in the heavy-top limit) Higgs exchange semi-inclusive e^e^ annihilation 
(SIA). Our results for the last quantity confirm the findings obtained in Ref. [9] by a different 
method, while the others are new. For brevity only a part of these results have been discussed here; 
for a full account the reader is referred to Ref. [18]. 

The numerical effect of the double logarithms beyond the third order on the spUtting functions 
in Mellin-A^ space is very small for a strong coupUng ots — 0.12 corresponding to a scale close to 
the Z-boson mass, amounting to less that 1 % for A'^ < 20. Here the contributions of terms beyond 
the fourth order are negligible. The corresponding corrections are larger, and receive noticeable 
contributions to order , for the coefficient functions. 

In (almost) all these cases these is no reason to believe that the N^LL corrections for £ > 2 
are small compared to the present results. This is hardly surprising, given previous experience 
with other end-point resummations, but indicates that more terms are required in order to achieve 
phenomenological relevance. We hope that the present results will provide useful information for 
the development of more sophisticated approaches in the future. 

As done at the fourth-order for non-singlet quark coefficient functions in Ref. [9] and the 
spacelike spUtting functions in Ref. [10], it is possible to extend the present results to all higher 
orders in the expansion in powers of (1— Jc), i.e., to all terms of the form (1— x)°ln^"~^~''o 
Consequently all large-x double logarithms in inclusive deep-inelastic scattering (DIS) and SIA are 
fixed by lower-order information, with the coefficient of the N^LL contributions determined by the 
N^LO fixed-order results. 

Due to the somewhat different structure of the D-dimensional phase-space integrations, the 
present approach is unfortunately not directly applicable (beyond the leading logarithms) to the 
Drell-Yan lepton-pair and Higgs-boson production in proton-proton collisions. Hence further, more 
refined tools are required in these cases to improve upon the physical-kernel constraints of Ref. [9] 
and to extend those results to the flavour-singlet contributions. 

The present approach can be extended, on the other hand, to high-energy (small-x) double log- 
arithms, if neither to all quantities nor to all powers a in the analogous small-x expansion in terms 
of x° In 2«-^-«o X terms. The NNLL resummation of the dominant x ~ ^ terms in the timelike splitting 
functions and the fragmentation functions in SIA has been presented in Ref. [33], the results for the 
terms for the corresponding DIS quantities have also been derived and will be presented else- 
where [34]. Clearly more research is required on both the threshold and high-energy logarithms. 
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